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Abstract. A theory of far-from-equilibrium transport in arrays of tunnel junctions 
is developed. We show that at low temperatures the energy relaxation ensuring tun- 
neling current can become a cascade two-stage process. First, charge carriers lose 
their energy to a bosonic environment via non-phonon energy exchange. The role of 
such an environment can be taken by electromagnetic fluctuations or dipole excita- 
tions (electron-hole pairs). The environment, in its turn, relaxes the energy to the 
thermostat by means of phonon irradiation. We derive the current-voltage charac- 
teristics for the arrays and demonstrate that opening the energy gap in the spectrum 
of the environmental excitations completely suppresses the tunneling current. The 
consequences of the cascade relaxation in various physical systems are discussed. 



1.1 Introduction 

Electronic transport in mesoscopic tunnel junctions is ensured by the en- 
ergy exchange between the tunneling charge carriers and energy reservoirs: 
since the electronic energy levels at the banks of the mesoscopic junctions 
are, in general, unequal, the tunneling is impossible unless there is a subsys- 
tem of excitations capable of accommodating this energy difference . In- 
tense studies of nano-structured and disordered systems including Josephson 
junctions ^, |o|, mesoscopic tunnel junctions Q and superconductors 
patterned superconducting films ||l^,|l4), highly disordered superconducting 
and semiconducting films [p^-pO[ reveal a prime importance of the out-of- 
equilihrium properties of an environment to which the tunneling charge carri- 
ers relax the energy. In particular, experiments on disordered superconducting 
films revealed that at extremely low temperatures the I-V characteristics ex- 
hibit highly nonlinear behavior indicating that the transport charge carriers 
are decoupled from the phonon thermostat. 
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Notably, the relaxation processes can occur via relaxation of energy from 
tunneling carriers to some other bosonic environment mediated the energy 
exchange between the tunneling current and phonon thermostat. The role of 
such bosonic mediator can be taken by ether the electromagnetic environ- 
ment [l|,P,^-p||2l|,|2^ by the neutral dipole excitations [the electron- hole (e-h) 
pairs] generated by the tunneling carriers , 24 . By consideration in the 



spirit of Feynman- Vernon influence functional it was shown [p5|-p7| that in 
the large Josephson junction arrays the low temperature transport is governed 
by the unbinding neutral dipole excitations and the successive tunneling of the 
resulting electric charges. In the course of tunneling the propagating charges 
generate the dipole environment controlling relaxation process. The formation 
of a gap in the dipole excitation spectrum impedes relaxation and causes the 
suppression of the tunneling current, i.e. the localization of charge carriers. 

The energy relaxation in mesoscopic tunnel junctions in the case where 
the energy exchange between the tunneling carriers and the electromagnetic 
and/or electron-hole reservoir, 1/Te_env, is comparable to the rate of the en- 
ergy loss to the phonon thermostat, l/rcnv^bath, was analyzed in [Q. In this 
article expanding on our kinetic approach developed in ||2^, we present 
a general approach to the description of strongly nonequilibrium processes 
where l/To_cnv ^ l/Tcnv^bath and show that the low-temperature energy re- 
laxation enabling the tunneling current occurs in two stages: (i) The energy 
relaxation from the tunneling charges to the intermediate bosonic modes, elec- 
tromagnetic or dipole (electron-hole) excitations, which we hereafter call the 
environment; and (ii) The energy transfer from the environment to the phonon 
thermostat, to which we will be further referring as to a bath. We demonstrate 
that at I/tc-ciw ^ l/Tcnv->bath, electronic transport is controlled by the first 
stage and is thus critically sensitive to the spectrum of the environmental 
modes. At the same time, the passing current drives the environment out of 
the equilibrium, and the environment spectrum and effective temperature may 
become bias-dependent themselves. We construct coupled kinetics equations 
for charge carriers and out-of-equilibrium bosonic environment and derive I-V 
characteristics for arrays of normal and superconducting junctions. 

The article is organized as follows. In Section 1.2 we construct a general 
theory of environment-mediated transport in a single tunnel junction and cal- 
culate the corresponding I-V characteristics. In Section 1.3, we extend our 
theory onto arrays of tunnel junctions; the content of these sections follows 
our recent publications [|9[0. In Section 1.4 we explore application of the 
cascade two-stage energy relaxation mediated transport to different physical 
systems and phenomena. In particular, we discuss nonequilibrium Coulomb 
anomaly, negative differential resistance and overheating in single junctions, 
properties of electron-hole environment in disordered metals, and the mecha- 
nism of tunneling transport in a superinsulating state re-deriving the results 
of [p5|-p7t from a different perspective. 
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1.2 A single junction 

We start with a tunnel junction between two bulk metallic electrodes biased 
by the external voltage V, see Fig. |l.l| a. A general formula for the tunneling 
current reads: 

I^e[f-f), (1.1) 

where 7^ (^) is the tunneling rate from the left (right) to the right (left), 
and, for a single junction, 

^ = i; L /i'Hi - ff)pHe - o , (1-2) 

where J^^'^' are the electronic distribution functions within the electrodes, 
P<(e) is the probability for the charge carrier to lose the energy e to the 
environment, and i?T is the bare tunnel resistance, representing the interaction 
of electrons with the bath. The backward scattering rate, F cx J^^, fe'^\l — 

fj:P)P^{e — e'). If an intermediate environment is absent and the relaxation is 
provided by the phonon bath, then P^(e) = (5(e) and Eq. ( p^ ) reproduces the 
conventional Ohm law. The quasiequilibrium situation where the distribution 
functions of the environmental modes Ni^ are Bose distributions parameterized 
by the equilibrium temperature was discussed in [^^,|3^. In a general, far from 
the equilibrium case, we find: 
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Fig. 1.1. (a) The effective circuit for a tunnel junction subject to bias V and 
with the environment having the impedance Z. (b)-(c) Diagrammatic expansion of 
to the first and the second orders in p, respectively. The solid lines represent 
propagation of electrons, the dashed lines denote the environment excitations. The 
vertex with the two electron lines and one dashed line carries a factor G-rpij^)/^, 
the two dashed-lines vertex corresponds to Gt p{u)) p{uj') / {ujuj' ) . 
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P<(e)==/ dtexp[J{t)+iet], (1.3) 

J —OO 

poo 1 

J{t) = 2 ^p[u)F{u), (1.4) 
Jo ^ 

F{uj) = [iV^e'"* + (1 + 7V^)e-''"* - B„] . 

Here exp[J(i)] is the nonequilibrium generalization of the Feynman- Vernon 
influence functional reflecting that tunneling electrons acquire random 
phases due to interactions with the environment, represented by a set of 
oscillators with the nonequilibrium distribution of modes, N^. The lat- 
ter is defined by the kinetic equation with the scattering integral describ- 
ing the energy exchange between environmental modes and tunneling elec- 
trons. Terms proportional to the and 1 -I- correspond to the ab- 
sorbed and emitted environmental excitations, respectively. The combina- 
tion = 1 + 2Ni^ is the kernel of the time-independent contribution to 
J describing the elastic interaction of the tunneling electron with the en- 
vironmental modes and having the structure of the Debye-Waller factor. 
In an equilibrium, N,^ reduces to the Bose-function and the functional 
recovers the result by Ref. The spectral probability of the electron- 

(electromagnetic) environment interaction is p{uj) — Re[Zt(a;)]/_RQ, where 
Zt = l/[iCuj + Z{uj)~^] is the total circuit impedance, Z is the environ- 
ment impedance, C is the junction capacitance, and Rq is the quantum re- 
sistance . Proceeding analogously to Ref. Q , one finds the spectral proba- 
bility corresponding to the electron-environment interaction within each elec- 
trode as Pn(<^) — 2 Im Un/ {Dnq^—iu)^, n — 1,2, and that for the interaction 

across the junction, pi2{u!) — — 2Im J^Uu/HDiq^ — iLu){D2q^ — iuj)], where 

£'i(2) are diffusion coefficients within respective electrodes, and f/i(2) are the 
dynamically screened Coulomb interactions within (across) the electrodes. 

The spectral probability p{uj) of the electron-environment interaction con- 
trolling the behavior of the I-V characteristics is central to our approach. The 
form of p{i-o), in its turn, is determined by the structure of the environmental 
excitation spectrum and, in general, depends on the external bias. The lat- 
ter dependence becomes essential in the array of highly transparent junctions 
where p{oj) is different for elastic and inelastic processes p],|^,p4| . In particular, 
for the e-h environment with constant U, one should cut off the (diverging) 
integral at q^^ = yJl/T^{Tce)D when calculating p(0). Here is the 
electron inelastic time and Toff is the (bias dependent) effective temperature 
of the environment defined below. This allows us to formulate a recipe: if in 
an equilibrium p = p{uj, T) then in an out-of-equilibrium state p = p{uj, Tcs). 
Of special importance are the effects of spectrum reconstruction accompany- 
ing the phase transformation in the system: it is the sensitivity of the shape 
of the I-V curves to the form of p = p{uj,T) that makes tunneling currents 
an irreplaceable spectroscopy tool. In particular, opening the gap T* in the 
environmental excitation spectrum suppresses p = p{uj) within the energy 
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interval < uj < T* giving rise to complete vanishing of the tunneling trans- 
port current at low temperatures T < T*. As we will discuss below, this is 
the microscopic mechanism behind the low-temperature suppression of elec- 
tronic transport in disordered arrays of tunnel junctions, the phenomenon of 
superinsulation. 

To close the set of formulas (p~l|)-( p^ ) one has to add the kinetic equa- 
tions (KE) for the boson distribution functions Nuj- To derive these KE we 
use a semi-phenomenological kinetic approach of [ p6[ and express the current 
of Eq. (LI) through the electronic distribution function as / = j^^ [dfiy' /dt]vi. 

Here ^1(2) is the density of states in the lead 1(2) and dfel^ /dt = Icoi, where 
Jcoi is the collision integral describing the evolution of the electronic distribu- 
tion function due to energy and/or momentum transfer processes. Expanding 
further with respect to p we obtain, in the zero order in N^^, the collision 
integral in a form I^^J = -J Wi2[f^l\l- f^?)- fS\l- f^l''mei-e2y2de2, 
where W12 = 1/viV2R-t is proportional to the bare probability for an electron 
to be transmitted from one lead to the other. In the first order 

^ = - / dLov^V2de2 {—\ W12X 
at J J 

Siei2 - coMIHn^ + 1)(1 - /if) - (1 - /if )A^./if ]+ (1-5) 

Sie^2 + c.)[/if A^Ul - /if) - (1 - ./if )(^^ + l)/if ] 

where ei2 = ei — £2 and i^uj is the density of environmental states ||37| . The 
structure of ^ifis identical to that of the electron-phonon scattering inte- 



gral in metals [ p6[ , where would stand for the phonon distribution func- 
tions. The quantity p/iLov^) is proportional to the probability of the electron- 
environment scattering. 

The collision integral dual to /i^j and describing the evolution of N^^ is 
derived analogously, and the resulting kinetic equation is: 

where A is the numerical factor of order of unity, is the electron-hole pairs 
distribution function. The scattering integral in Eq. (^^) is also identical by 
its structure to the phonon-electron scattering integral in metals [|6| . For the 
electron- hole environment [i — 1,2 labels the electrode in which the pair is 

located), one has nL'^ = (l/w) /i*^(l - /ii^), where e± = e±a;/2; this agrees 
with the results of Ref. p7| where the nonequilibrium bosonic distribution 
function is equivalent to our 1 -|- 2nui. If electrons and holes belong to different 
electrodes, 

n(;^-) = (2a.)-i //«a,^-(l-/(f), (1.7) 
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where is the Pauh-matrix. From Eq. (1.6) one estimates the rate of 
the energy exchange between the environment and the tunnehng electrons 
as: l/re_env = p{'-^) / iy^^R^) . Now one has to compare l/ro_cnv with the 
rate of the interaction of the environment modes with the (phonon) bath, 
l/7'cnv^bath('^)- For the electron-hole environment, l/Tenv->bath('^) is deter- 
mined from Eq. ( |l.6| ) to which the electron-phonon scattering integral is added. 
If Tenv->bath ^ Tc-onv, the two-stagc energy relaxation takes place and the 
characteristic energy transfer from tunneling current is a; ^ max{Tc,F}, 
where is the electronic temperature in the leads. The electromagnetic 
environment mediates the two-stage relaxation in the case where Ohmic 
losses occur in a LC superconducting line and are small To take a 

typical example, in aluminum mesoscopic samples To_env = 10~* sec and 
Tenv->bath = 10~^ scc [|2| , SO the Conditions for the two-stage energy relax- 
ation are realized. Then the distribution functions, iV^, deviate significantly 
from the Bose-distribution with the temperature of the phonon bath. They 
should be determined from the condition that the collision integral of the 
environmental modes with the e-h pairs accompanying the current flow be- 
comes zero, and Eq. ( |l.6| ) yields iV^ = n^'^\ If To ^ V, then iV^^ can be 
approximated by the Bosc-function with some effective temperature T^s at 
uj < V — Teff and N^^ = at w > Tcff (the emission of the excitations with 
the energy larger than V is forbidden), and 

Tcff EE lim LoN^ = O.bV coth(F/2T) . (1.8) 

a;— ^0 

This result shows that the system with the environment well isolated from 
the bath cannot be cooled below Tcff- Note that the coth-expression for Toff 
obtained in the first approximation in p. In a general case. Tog depends on p. 

Equations ( |l.l| )-( pr6| ) give the full description of the kinetics of the tunnel 
junction in a nonequilibrium environment. To derive the I-V characteristics we 



find iV^ = n-u"^' and plug it into Eqs. (1.1)- (1.4). Introducing the parameters 



g^^ — p{Q) and A, the characteristic frequency of the p(w) decay [for the 
Ohmic model [Q, p = / {\ + (w/yl)^} and A/g is of the order of the 
charging energy of the tunnel junction], we find: 

V A 



in the interval T < < yl, where Tog ~ V. Note that I{V) given by Eq. (1.9) 



differs from the power law dependence obtained in 31 1 for — Tog = 0. This 
shows that tuning the environment one can control the /(y)-characteristics 
of the tunnel junction (the gating effect). At high voltages, V A, one finds 

I{V)~[V ~ A^)/R^, (1.10) 

poo 

Aoo = ^ J'(0) = 2 / dLop^[l + 7V(°"*) - iV(-)] ~ Z\W ln(yl/ min{r„ Te„v}), 

(1.11) 
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Fig. 1.2. (a) The single electron two- islands' circuit, (b)-(e) Diagrams describing the 
forward inelastic cotunneling rate. The "up" arrows stand for the e-h pairs excited 
during the cotunneling and the "down" arrows correspond to the recombination of 
the e-h pairs. The vertices shown by boxes are proportional to the probability of an 
elemental e-h pair excitation, p{uj)/uj. 



where A 



(0) 



(in) 



Ao^[N'^°''^^ = iVC")] - Ajg, since at V" » yl, 7Vi°"*^ c^; Ajuj > 



1.3 Arrays of tunnel junctions 

Extending Eq. (|l.2[) onto an array comprised of N junctions one finds 



/ ^ 



Rc 



/ dede'h{e)[l~h{e'me^e'), 



(1.12) 



where 



dtexp{iEt)\ / duj 

-OO I Jo 



X 



r(in) iut 



(out)x -icjt 



j<N~l 



(1.13) 
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Here S ^ Ec '^^^^^N^/{N - 1)!, and = e^/2C is the Coulomb charging 
energy of a single junction (C is a single junction capacitance) and for the 
Cooper pair transport e 2e. Equations ( |l.l2[ ) and (1.13) were derived in 
a first order in tunneling Hamiltonian. Shown in Fig. |l.2| is a diagrammatic 
representation of Eq. ( p_.13| ) for = 3. 

A generalization of the results obtained for a single junction including the 
structure of the collision integral and the concept of the effective temperature 



Eq. (1.8), onto large arrays is straightforward. As long as temperatures are 
not extremely low [ p3[ , the charge transfer in large arrays is dominated by 
the inelastic cotunneling and the cascade energy relaxation. The tunneling 
carriers generate e-h pairs [^,|2^ serving as an environment exchanging the 
energy with the tunneling current and then slowly losing it to the bath. 



1.4 Cascade two-stage relaxation: general applications 

In the preceding sections we have formulated a general approach to descrip- 
tion of electronic transport in mesoscopic tunnel junctions mediated by the 
energy exchange between the charge carriers and the environment, which, 
in its turn, relaxes the energy to thermostat. In what follows we will show 
that there is a rich variety of the seemingly disparate kinetic phenomena 
in tunnel junction arrays that allow for a natural and transparent descrip- 
tion within the unique framework of this hierarchical or cascade relaxation 
concept. We illustrate the power of our approach applying it to phenomena 
of overheating |38, 3^, the Coulomb anomaly the transport in the 



granular systems governed by the cotunneling processes |31|, |40|] and medi- 
ated by the electron- hole environment. The concept of the cascade two-stage 
energy relaxation appears to be of the fundamental importance to revealing 
the microscopic mechanism of the insulator-to-superinsulator transition and 
the nature of low temperature transport in the superinsulating state |p5|-p7|| . 
We show, finally, that the concept of the cascade energy relaxation resolves 
the long standing puzzle of temperature-independent pre-exponential factor in 
variable range hopping conductivity |l"^-|r^,|4l|-^. We demonstrate that the 
underlying physical mechanism behind all the above is the imbalance between 
the intense energy exchange of the charge carriers with the (nonequilibrium) 
environment and the comparatively weak coupling of the environment to the 
phonon bath. 



1.4.1 Nonequilibrium Coulomb anomaly 

Consider the quasi- 2D disordered metallic film attached at the edges to two 
electrodes, see Fig. L3(a). The Coulomb anomaly arises due to electromag- 
netic fluctuations associated with the electron-electron interaction in the wire 
(electron-hole environment excitation) and manifest itself as the suppression 
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(b) 



^(1) 
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Fig. 1.3. (a) The quasi-2D disordered metallic film attached at the edges to two 
electrodes. The tunneling density of states (TDOS) can be determined using the 
transport measurement done with the help of the Scanning Tunnel Microscope 
(STM) tip. (b) Equivalent circuit that describes transport between the STM-tip 
and metal (in nonequilibrium state). 



of the local tunnel density of states (TDOS) at small energies [Q. Experimen- 
tally the density of states can be determined by means of tunneling transport 
measurements with the scanning tunnel microscope (STM) tip serving as one 
of the electrodes, see Fig. 1.5(a) and Refs. |ll|,|2^]. The corresponding equiv- 
alent circuit is shown in Fig.^]j(b). The resistor characterizes the contact 
between the tip and the quasi-2D metallic film. This is the same resistor that 
appears in Fig. |it| . 

In the absence of the current, i.e. in an equilibrium, the distribution func- 
tion of electrons at any point of the film is the Fermi function. Then the 
current voltage characteristics of the junction, /(Vtip), where Vtip is the po- 
tential of the tip, can be found following the recipes of Refs. we have 
to use Eqs. ( |l.l| )-( pr2| ) with the equilibrium environment distribution function 
and take p{uj) = 2Im U /{Dq^ — ioj)"^, where U is the dynamically-screened 
Coulomb interaction in the metal and D is the diffusion coefficient. Then from 
the differential conductance, dl/dVtip, we get the local density of states (the 
same as in Ref. [Q) and find all the standard Coulomb anomaly features, see 
Refs. 

As the current starts to flow through the metallic film, then the distri- 
bution function of electrons, /(e, r), at low temperatures (where the phonon 
bath is frozen out) becomes nonequilibrium and should be found from the 
kinetic equation ||ll|,|2^. A/ = with the edge conditions relating electron 
distribution function in the metal with the electron distribution functions at 
the reservoirs. The solution of the kinetic equation at the center of the metallic 
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film can be approximated as follows: 



(1.14) 



where V is the electrical poten tial of the left electrode that pushes the cur- 
rent through the metal, see Fig. 1.3(a). Here it was assumed that the electron 
inelastic scattering length, ^in, is larger than the separation between the elec- 
trodes attached to the metallic film. The electron-hole environmental modes 
in the metallic film, which are responsible for the Coulomb anomaly, become 
now nonequilibrium. Since the phonon bath is frozen, then Tcnv-bath 3> ro_env, 



so the stationary distribution function A'^j is to be found from Eq. (l.C). Then 
we obtain, 



(2^) 



'J de/e+(l-/.. 



(1.15) 



We calculate the tunneling rate between for electrons passing from the tip to 
the metallic film using Eq. (1.2): 



7^ = ^ 



/fp)(l-/,OP<(e-e'), 



(1.16) 



where /i'''^'' — /i?(e — Vtip) and is given in Eq. ( 1.14 ). Note that when eval- 
uating P^, we used the nonequilibrium environmental bosonic distribution 
function N^j as given by Eq. ( 1.15 ). Equations (1_^) and ( 1.16 ) yield the cur- 
rent, flowing through the tip, /(Vtip), from which one determines the TDOS 
as: 



dI{Vtip) 



dVtip 



(1.17) 



One finds after some algebra that the TDOS acquires the 'superstructure' 
on the energy scale of the order of Tcs{V): two dips develop at the energies 
corresponding to the positions of the Fermi level in the tip and in the film, 
where the electronic distribution function experiences discontinuous jumps 
(i.e. at e = and e = y) |ll|,||]. 



1.4.2 Negative differential resistance and overheating 

If the rate of the energy supply from the external bias to the charge carriers 
exceeds the rate of energy losses to the environment, the phenomenon of 
overheating takes place and the energy distribution function of the current 
carriers noticeably deviates from the equilibrium distribution function [^ , |39[| . 
One of the characteristic manifestations of the overheating effect is the onset 
of the "falling" region of the I-V curve where the differential conductivity 
G = dl/dV < 0. The corresponding I{V) characteristics is referred to as that 
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of the S'-type if the current is the muhi-valued function of the voltage, and 
the I-V curve of the A'^-type corresponding to the case where the current is 
nonmonotonic but still remains a single valued function of the voltage. The 
phenomenon of overheating has been a subject of the incremental interest 
and extensive studies during the decades, see the detailed Volkov and Kogan 



review |38 , and the impressive progress in understanding of the underlying 
mechanisms was achieved. Recently, for instance, the ideas of [ p8[ were applied 
to description of nonlinear low temperature current- voltage characteristics in 
disordered granular metal [Q. The concept of the two-stage relaxation of 
the present work enables us to construct a general approach to overheating 
taking into the consideration also non-phonon mechanisms of relaxation. A 
general scheme is as follows. Let us consider the far-from-equilibrium tunneling 
transport mediated by an environment. Let the environmental excitations 
emitted by tunneling carriers be "hot" ("up" lines in Figs. p^|l.2[ ) whereas the 
environmental excitations absorbed by tunneling carriers are "cold" ( "down" 
lines in Figs. |l.l[ - |l.2| ). This implies that there exist two distinct energy scales 
characterizing the frequency distributions, Tvi,'"'' and ivi°"*\ corresponding to 
emitted and absorbed environment excitations, respectively: 

T^^r^ = lim '^A^^°"*\ (1.18) 

CO — 

ri'''^ = lim wM'"). (1.19) 

These two temperature scales are, in general, different provided the effec- 
tive inelastic length, L^, describing thermalization of the excitations with the 
phonon bath (thermostat) is of the same order or smaller than the hot region 
where tunneling electrons strongly interact with the environment. Then the 
absorbed environmental excitations come from the outside of the hot region 
and have the temperature of the bath, T^g"^ w Jbath- The emitted excitations 
are hot, and their effective temperature is the function of the applied voltage 
(current). The temperature of the hot excitations is to be found from the heat 



balance equation following the recipes of Refs. |3q, 39 . Namely, it is obtained 



by integrating the product of the kinetic equation for electron distribution 



function [see, e.g., Eq. (1.5)] and the electron energy e, over the volume and 
energy: 

_ ,'7.(111) \ .„(0Ut)s' • > 

'c— cnvV-^cff / 'c— cnvV-^cff / 

where ro_onv(T') is the temperature dependent electron-environment inelastic 
scattering rate, £{T) is of the order of total the energy of electrons having the 

temperature T in the volume involved and P F^/i?(T^'™*'') is the (Joule) 
heat produced by the passing current. The I-V characteristics comes out as 



solutions to Eq. (1.20) with the specific relations between £ and the applied 



voltage and current. The resulting I-V curves may come up as the so called S- 
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or TV-type of the I-V characteristics depending on the specific temperature 
dependencies R(T^^^^^), see Ref. 



. Making use of Eq. (1.5) we can esti- 
mate Te_env(^r). It is proportioual to the coefficient standing in the scattering 



integral by /i^^ 



<5(612 - ^)[{N^ + 1)(1 - /if) + W./if ] + 
5{e,2+Lo)[N^{l-I^f) + {N^ + l)fif] 



(1.21) 



The equation (1.21) is derived for a single tunnel junctions. To describe 
the overheating effects in arrays of many junctions one needs an appropri- 
ate generaUzation of a single junction approach. An important example of 
such a multi-junction system is an array of metallic granules in the insulating 
state. In granular metals the role of environment is taken by the electron-hole 
pairs. In this case To_onv is electron inelastic relaxation time, T^^'^ = T is 
the bath temperature. The energy scale £ can be roughly approximated by 
the energy of the quasiparticles in the bulk of the sample of the volume J7 
as S{T) ^ ir^i/SlT^ /6. Employing the standard theory of the tunneling con- 
ductivity in the granular materials [^4| one then can arrive at the S'-type 
/-1^-characteristics for granular metals. 



1.4.3 Electron- hole environment in a disordered metal 

To gain an insight into the behavior of the electron-hole environment in meso- 
scopic tunnel junctions and reveal the role of Coulomb interactions, which 
are instrumental to formation of the environment properties, let us employ 
the Finkel'stein's theory of Fermi-liquid fluctuations in a disordered granular 
metal . The theory treats the Coulomb potential as the contact interaction. 
This offers a pretty good description of disordered metals, where Coulomb 
potential is well screened to i5-function. Of course, in the close vicinity of 
the disorder-driven metal-insulator (or superconductor-insulator) transition 
where the real Coulomb fields come into play, the orthodox Finkel'stein's the- 
ory is not valid, and one needs the approach capable to accommodate the 
long-range Coulomb effects. 

We adopt the Keldysh representation which is most adequate for discussing 
dynamic effects. The action describing Fermi-liquid fluctuations acquires the 
form: 

iSp — — tr 

(1.22) 



iTp 



{iPi)txip2)tx + 2ii^tr[^''(Tx 0]} 
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Here Pi(2) is the operator of the Fermion fluctuations density in the classical 
(quantum) sector, 



Pi 



2ttv 



tr(o-xQt.t') + TT" 



P2 



27riy 
1 + F„ 



tr(a-x(5t,t') + ^ 
Ztt 



are the conjugated fields, Fp — Fp/{1 + Fp), Fp 



(1.23) 

is the contact ampli- 



tude modeling the screened Coulomb interaction in the singlet channel. The 
Fourier transform of the typical retarded fiuctuation propagator describing 
fluctuations has the form: 



2?" = 



1 



Dq^ - iLj{l + Fp) ' 



(1.24) 



where uj are the Matsubara frequencies. The diffusion poles of the propagators 
thus possess the structure, Dq'^ — icu, Dq^ — jw(l + Fp). One sees that this 
description applies and the spectrum of the excitations is stable only as long as 
the single granule mean level spacing remains the smallest energy parameter, 
i.e. S < Dq^, uj, ujFp. This implies the development of the excitations spectrum 
instability and opening the energy gap at temperatures 



T <T* = — 

J- n 



(1.25) 



Formation of the energy gap in the spectrum of the environmental excitations 
is a phenomenon that tremendously influences the tunneling electronic trans- 
port. One sees straightforwardly that as the energy gap, T*, in the spectrum 
of electron- hole excitations appears at T < T*, the spectral probability p{uj) 
for the electron-environment interaction vanishes in the interval < w < T*. 



Then Eq. (1.12) yields the suppression of the tunneling current to the prac- 
tically zero magnitude. The physical significance and meaning of this result 
is merely that as the gap in the excitation spectrum opens, the environment 
ceases to efficiently mediate the energy relaxation from the tunneling carriers, 
impeding thus the tunneling current. 

Interestingly, the similar temperature that marked vanishing of the con- 
ductivity by the e-e interactions in the absence of coupling of electronic system 
to phonons in the model of disordered quantum wire was found in Refs. [ ^|52| 
by reducing the electron conductivity to the Anderson model on the Bethe 
lattice. This suppression of conductivity was interpreted as the Anderson lo- 
calization in Fock space. We reiterate here that models based on the contact 
e-e interactions hold only as long as the long-range Coulomb effects are effec- 
tively screened. As a result its applicability, for example, to the description 
of physical phenomena in the vicinity of metal-insulator or superconductor- 
insulator transition needs special justification. In the next section we propose 
an approach capable to explicitly account for the long-range Coulomb in- 
teractions in the critical region of disorder-driven superconductor-insulator 
transition. 
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1.4.4 Superinsulating behavior 

Let us consider charge transfer in a two-dimensional array of superconducting 
tunnel junctions (or, equivalently, Josephson junction network) in the insulat- 
ing state, i.e. under the conditions that Ej < Ec, where Ej is the Josephson 
coupling energy, and Ec is the charging energy related to the capacitance C 
between the two adjacent granules (or the capacitance of a single Josephson 
junction). We focus on the limit C ^ Cq, where Co is the capacitance of a sin- 
gle junction to the ground. The electric properties of the array are quantified 
by the screening length A = ay^C/Ca, where a is the size of the elemen- 
tal unit of the Josephson junction network. At distances R < X, the electric 
charges interact according the logarithmic law, the energy of the interaction 
being cx ln(i?/A), at larger distances Coulomb interaction is exponentially 
screened. If the linear dimension of an array, L does not exceed A, the charge 
plasma, comprising the environment in the superconducting array, experiences 
the charge binding-unbinding Berezinskii-Kosterlitz-Thouless (BKT) transi- 
tion [|3[-|8| at T = Tbkt ^ Ec |^,|5|. In the low temperature BKT phase 
positive and negative charges are bound into the dipole pairs, whereas above 
the transition the charges form a free plasma. Binding positive and negative 
charges into pairs implies that aX T :^ E^ the energy gap ~ Ec opens in the 
spectrum of unbound charges. As we have discussed in the preceding section, 
opening the energy gap in the environmental excitation spectrum results in 
vanishing p{ijj) in the interval < u! < Ec. 

The low temperature transport in the large JJA arrays was discussed 
in [?,^, ^ within the framework of the Feynman- Vernon influence func- 
tional technique. The charge transfer is realized by the dipole excitations that 
unbind and further tunnel under the applied bias. In the course of tunneling 
the propagating charges generate and leave behind the strings of the dipole (or 
electron-hole) excitations that accommodate the charges energy difference at 



different superconducting islands. As shown in 25-27]] the system of dipoles 
excited at Josephson junctions can be mapped onto the ensemble of quan- 
tum rotators. The subsequent averaging of the current with respect to all the 
states of the rotator ensemble is equivalent to averaging over the thermody- 
namic bath, the role of which is thus taken by the dipole environment. Thus 
the relaxation process ensuring the tunneling current is governed by the dipole 
environment and the subsequent phonon emission by these dipole excitations. 

Now we discuss the same process from the viewpoint of the nonequilib- 
rium approach developed in the present work. Following further the line of 



reasoning used in preceeding sections and employing again Eq. (1.12), we ar- 
rive at the conclusion that at T < Tbkt, the charge transfer in the array of the 
superconducting tunnel junctions is suppressed. Remembering now, that at 
low temperatures tunneling current in an array of superconducting junctions 
is governed by the cotunneling processes |2^,|2^ described in Section 1.3. An- 
alyzing the contribution from higher orders into the cotunneling process, one 
finds that the current suppression holds in all orders. 
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On a qualitative level the effect of the gap opening in the charge plasma ex- 
citation spectrum on the charge transfer in the superconducting tunneling ar- 
ray can be described as follows. Starting with Eq. (1.1) and using Eqs. ( 1.12| )- 
( 1.13| ), one can find 

/ oc exp(-£;/H/) , (1.26) 

where E is the characteristic energy for the charge transfer between the gran- 
ules. In the system with the linear size L < X, E = Ecln{L/a), where 
Ec = 2e^/C for the Cooper pair, and E^ = e^/2C for the quasiparticle. 
The quantity W is the energy scale associated with the rate of the energy 
exchange between the tunneling charges and the environment. While the rig- 
orous derivation for W is not available at this point, the estimates suggest the 
interpolation formula 

W ~ . (1.27) 

exp(£;e/T)-l ^ 

The meaning of this formula is that the relevant energy scale characterizing 
the tunneling rate is the energy gap that enters with the weight equal to 
the Bosc-kind filling factor describing the probability of exciting the unbound 



charges. Well above the charge BKT transition Eq. (1.27) gives W = T (this 
reflects the fact that the number of the unbound charges is determined by the 
equipartition theorem and is proportional to T /Ec). This yields 

/ cx exp[-£;cln(L/a)/r] , T>T^^^. (1.28) 

One has to bear in mind that this formula holds only at temperatures not 
too high above the charge BKT transition, where Coulomb interactions are 
not completely screened. At T ^ TeKT, E ^ Ec ln(A/a) and as the screening 
length becomes of the order of a, E Ec- In this temperature region the sys- 
tem exhibits the 'bad metal' behavior. Notably, Eq. ( 1.2^ ) looks like a formula 



for the thermally activated current, however one has to remember that the 
physical mechanism behind the considered charge transfer is quantum me- 
chanical tunneling process which can take place only if the mechanisms for 
the energy relaxation are switched on. 

At low temperatures, T <^ Tbkt — Ec, the characteristic energy from 



Eq. ( 1.27 ) isW — Ec exp(— E'c/T). The unbound charges that have to mediate 
the energy relaxation from the tunneling carriers are in an exponentially short 
supply, and one finds 

/ocexp[-ln(L/a)exp(£;c/r)], T < Tbkt , (1-29) 



reproducing the results of |25|-|27[| for the Cooper pair conductivity. The ad- 
ditional suppression of the electronic transport in Josephson junction arrays 
as compared to the activation regime was indeed experimentally observed in 
Refs. 1591,^. 
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The outlined picture of the tunnehng current in the superconducting junc- 
tion array applies to thin superconducting films close to superconductor- 
insulator transition (SIT) p6|-p^. Indeed, in the close vicinity of the SIT, 
the disorder-induced spatial modulations of the superconducting gap give rise 
to the electronic phase separation in a form of the texture of weakly cou- 
pled superconducting islands with the characteristic spatial scale of the 
order of several superconducting coherence lengths ^. The consequences of 
this phase separation are two-fold. First, the fine balance between disorder 
and Coulomb forces near the SIT results, on one hand, in the effective re- 
duction of the disorder strength, and, on the other hand, in the emergence 
of the strong non-screened Coulomb field in the spaces between the islands. 
This field can be argued to reconstruct the island texture into a nearly reg- 
ular array of weakly coupled superconducting islands, i.e. into an array of 
mesoscopic superconducting junctions. Note that the regular lattice of the 
superconducting islands can be also induced by nonlocal elastic fields due to 
coupling between the film and the substrate . Second, in the critical region 
near the SIT this island array is on the verge of the percolation-like transition 
between the superconductor and insulator (the cluster of the islands touching 
each other and traversing the sample means the emergence of global supercon- 
ductivity). Since in such two- phase systems the dielectric constant diverges on 
approach the percolation transition |63|-p7| , the thin films in the critical 
region of the SIT develop a huge global dielectric permeability. This means 
that on the distances L < ed, where d is the film thickness, all the electric 
fields, emerging due to fluctuations in the local charge of the 2D environment, 
remain trapped within the film, and, therefore, the charge environment is 
nothing but the 2D neutral plasma, where charges interact with each other 
according to the logarithmic law. Therefore the 2D charge environment expe- 
riences the binding-unbinding charge BKT transition at T = Tbkt, and the 
energy gap opens in the environmental excitations spectrum at T < Tbkt, 
giving rise to suppression of both Cooper pair- and normal quasiparticle cur- 
rents, since both tunneling processes - irrespectively to whether it is transfer 
of Cooper pairs or quasiparticles across the junction - can occur only if the 
energy exchange with an environment is possible. Thus opening the gap in 
an environmental spectrum due to long-range Coulomb effects and the re- 
sulting suppression of the tunneling current offers a microscopic mechanism 
behind the insulator-to-superinsulator transition and the conductivity in the 
superinsulating state | 

Note in this connection, that in the recent work [|68[ the dc conductivity of 
an array of Josephson junctions in the insulating state is discussed in terms of 
transport of Cooper pairs through the narrow band formed by the Josephson 
coupling and distorted by weak disorder. This model for the Cooper pair 
transfer is in a sense complimentary to that of the present work where we 
discuss the opposite limit of sequential tunneling between the granules with 
the essentially different tunneling levels. 



25-27 
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1.4.5 Variable Range Hopping conductivity 

The notion of the cascade two-stage relaxation is a key to resolving the contro- 
versy of the variable range hopping (VRH) conductivity. Many experimental 
studies of the hopping conductivity in doped semiconductors p2|-44|, in ar- 
rays of quantum dots and in disordered superconducting films [16- ^^ , 
revealed the temperature independent pre-exponential factor, indicating the 
non-phonon mechanism of relaxation. What more, it had the universal values 
of the integer multiples of e^/fe. According to the early ideas by Fleishman, 
Licciardello, and Anderson [|69|, the universal pre-exponential factor may ev- 
idence that the energy relaxation is due to electron-electron (e-e) rather than 
the electron-phonon interactions. On the other hand, according to Refs. 
the e-e relaxation alone (i.e. in the absence of the coupling to phonons, which 
are in any case present in a real physical system) cannot ensure a finite con- 
ductivity below the certain localization temperature. Namely, in the absence 
of phonons the system of interacting electrons subject to disorder undergoes 
a localization transition ('many-body localization') in Fock space. The con- 
cept of the cascade energy relaxation developed in the present work resolves 
this controversy. The first stage is the fast energy exchange between the hop- 
ping electrons and the electron-hole environment (electrons do not see sparse 
phonons). Since this process involves only electron-electron interactions, the 
resulting pre-exponential factor does not depend on temperature. The second- 
stage process is the transfer of the energy from the electron-hole environment 
to the phonon thermostat. This process which is of course necessary to ensure 
the current, does not influence the pre-exponential factor. 



1.5 Conclusions 

In conclusion, we have developed a theory of the far from the equilibrium 
tunneling transport in arrays of tunnel junctions in the limit l/rc_env ^ 
l/'T'onv-bath- Wc havc demonstrated that the energy relaxation ensuring the 
low-temperature tunneling current occurs as a cascade two-stage process: the 
tunneling charges lose their energy to an intermediate agent, environment, 
and the latter relaxes the energy further to the thermostat. We have de- 
rived the tunneling I-V characteristics and shown that they are highly sensi- 
tive to the structure of the spectrum of the excitations of the environmental 
modes. In particular, opening the energy gap in the excitation spectrum sup- 
presses the tunneling current. As an important example we have considered a 
two-dimensional array of the superconducting junctions where the two-stage 
relaxation occurs via the two dimensional Coulomb plasma of positive and 
negative charges with logarithmic interaction, which experiences the charge 
binding-unbinding BKT transition at T = Tbkt- We argued that the gap in 
the plasma excitation spectrum that emerges gives rise to sup- 

pression of both, tunneling Cooper pair- and quasiparticle currents, thus of- 
fering the possible microscopic mechanism for the insulator-to-superinsulator 
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transition and the low temperature transport in the superinsulating state. We 
considered applications of our general approach to several physical systems 
and low temperature transport phenomena, including Coulomb anomaly in 
2D disordered metals, overheating effects, electron-hole environment in dis- 
ordered metals, and the origin of temperature independent pre-exponential 
factor in hopping conductivity. 
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